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ABSTRACT 

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic 
Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles 
approximately centered at /=30° and 1=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a 
widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, 
for a good fraction of them where the distance could be estimated. The compact sources (which we will call cores' in the following) are found for 
the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a 
core seems to appear when a threshold around A v ~ 1 is exceeded for the regions in the /=59° field; a A v value between 5 and 10 is found for the 
/=30° field, likely due to the relatively higher distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into 
filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a 
few for many sources, we find core surface densities between 0.03 and 0.5 g cirT 2 . Our results are in good agreement with recent MHD numerical 
simulations of filaments forming from large-scale converging flows. 
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1. Introduction 

From the diffuse cirrus to the molecular clouds, onto the forma- 
tion and death of stars, the Galactic Plane is the set where all the 
phases of the Galaxy life-cycle can be studied in context. Dust, 
best observed in the infrared and in the submillimeter, cycles 
through all these phases and is, as such, a privileged tracer for 
the Galactic ecology . IRAS ( |Neugebauer et al.|[l9iS) l and COBE 



(Mather et al.[|1990) were of tremendous importance in boosting 



the research in Galactic star formation and interstellar medium to 
the prominent positions they have today. As remote as they may 
now seem, however, these missions are only some 20 years away. 
Since then, a continuing explosion of Galactic Plane surveys, 
both in the mid-infrared at A < 70/itn (Omont et al. 2003 Price 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



et al.|2001 1 [Benjamin et al.|2003||Carey et a l. 2009) and in the 
submillimeter at A > 800jt/m (Schuller et al. 2009 Rosolowsky 
et al. 2009 ), are assembling a picture where the Galactic Plane 
has become accessible at sub-30" resolution over three decades 
of wavelength. The exception is the critical interval between 70 
and 500/mi where the bulk of the cold dust in the Galaxy emits 
and reaches the peak of its Spectral Energy Distribution (SED). 
The Hi-GAL Key-Project (Herschel infrared Galactic Plane sur- 
vey) will fill this gap. 

Hi-GAL is the key-project (KP) of the Herschel satellite 
(Pilbratt et al.| 2010} that will use 343 hours observing time 
to carry out a 5-band photometric imaging survey at 70, 160, 
250, 350, and 500/im of a | b |< l°-wide strip of the Milky 
Way Galactic Plane in the longitude range -60° < / < 60°. Hi- 
GAL is going to be the keystone in the multiwavelength Milky 
Way, opening up unprecedented opportunities with a promise of 
breakthroughs in several fields of Galactic astronomy. A full de- 
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Fig. 1. Three-color image (blue 70fim, 
350/mi) of the 2°x2° field around Z=30° 



green 160/im, red 



Fig. 2. Three-color image (blue 70/im, green 160/im, red 
350//m) of the field around Z=59° 



scription of the survey and its science goals are given elsewhere 
(Moli nari et al.[[2010b[ ). This contribution presents the first Hi- 
GAL data obtained in the Herschel science demonstration phase 
(SDP) and describes a few of the main early results that will be 
detailed in other contributions in this volume. 



the sensitivity predictions for point source sensitivity from the 
HSpot time estimator for all bands except at 70/im where the 
predicted limit is effectively reached, confirming that the noise 
in our maps is dominated by the cirrus confusion at all wave- 
lengths. A more detailed quantitative analysis is presented by 
|Martin et al.| ( [20T0| >. 



2. Observations and data reduction 



The H erschel PACS ( Poglits ch~eTaT| [20T0| and SPIRE (Griffin 
|et al.| |2010| l imaging cameras were used in parallel mode at 
60"/s satellite scanning speed to obtain simultaneous 5-band 
coverage of two 2° x 2° fields approximately centered at [I, 
b]=[30°, 0°] and [59°, 0°]. The detailed description of the ob- 
servation settings and scanning strategy adopted is given in 
Moli nari et al.| (|2010b). Data reduction from archival data to 
Level 1 stage was carried out using the Herschel Interactive 
Processing Environment (HIPE, Ott 2010) using, however, cus- 
tom reduction scripts that considerably departed from standard 
processing for PACS (Poglitsch et al. , 2010) and, to a lesser ex- 
tent, for SPIRE ( [Griffin et al.||2010| ). Level 1 Time Ordered Data 
(TODs) were exported from HIPE into FITS files. Further pro- 
cessing including the map generation was carried out using ded- 
icated IDL and FORTRAN codes. Saturation conditions were 
reached for all detectors only in SPIRE 250/im and 350/im im- 
ages in correspondence with the 3 brightest peaks in the /=30° 
field. The prescribed flux correction factors for PACS (Poglitsch 
etaL][20T0l ) and SPIRE ( jSwinyard efaL] [20T0] ) were applied to 



the maps since their photometric calibration was carried out us- 
ing the default calibration tree in HIPE. A detailed description 
of the entire data processing chain, including the presentation of 
the maps obtained in the five bands for the two observed fields, 
can be found in Traficante et al. ( 2010[ ). In the present letter we 
present in Figs. [T] and [2] the three-color images obtained using 
the 70, 160, and 350/im data (7=30° and 1=59°, respectively). 

These amazing maps convey the immediate impression of 
extended filamentary structures dominating the emission on all 
spatial scales. Measurements of the standard deviation of the 
signal at all wavelengths in the lowest brightness regions of 
the 1-59° field yield average values a factor two higher than 



3. Results and science highlights 

3. 1 . From IRDCs to mini-starburst and their impact on the 
ISM 

HerscheTs ability to observe such large areas with unprece- 
dented wavelength coverage and extraordinary signal dynamical 
range allows us to image simultaneously progenitors clouds for 
massive protoclusters to entire clusters of Young Stellar Objects 
(YSOs) in acitve star forming regions, while also measuring the 
effect of their strong stellar winds and powerful outflows on the 
surrounding medium. 

Infrared Dark Clouds (IRDCs) have received considerable 
attention in recent years (e.g. Rathborne et al. 2006 and Peretto 
& Fuller 2009 ) as potential sites for precursors of cluster forming 



sites. Found in silhouette against the bright mid-IR background, 
they shine in emission with Herschel. Peretto et al. (2010) shows 
how temperature effectively decreases from ambient values (20- 
30K) down to T=8-15K toward the densest (~ 10 23 cirT 2 ) peaks 
of these objects, resolving further temperature substructures that 
can be proxies for subsequent fragmentation. 

At the other end of the massive star formation timeline, we 
find W43, visible in the left portion of Fig. [T] as an outstand- 
ing case of Galactic mini-starburst. Detailed SED construction 
and luminosity estimates allow us to assess the very early evo- 
lutionary stage of the most luminous and massive YSOs in the 
region. It is remarkable how the same images show a promi- 
nent ridge extending southward which encompasses a 70pc-wide 
large cavity excavated by the W43 cluster and which possibly 
triggers further star formation (Bal ly etaL] [2010). Triggered star 
formation in less extreme environments can also be studied in 
statistically significant fashion modeling the SED of the sources 
found in correspondence of the multitude of Hn-driven bubble- 
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like structures found in the images, as shown for the bubble N49 
by |Zavagno et al.| ( [20T0| . 

Feedbacks from massive star formation, together with the in- 
tricate relationship between the interstellar radiation field and 
molecular clouds, are at the origin of the observed complexity 
of the ISM emission structure, where temperature ranges from 
~10K of pre-stellar cores to the ~40-50K of the photodissocia- 
tion regions ( |Bernard et al. 2010) >. 



3.2. Census of compact sources 

The extraction of compact sources is quite a challenging task in 
these fields, which we faced using a novel approach based on the 
study of the multidirectional second derivatives in the image to 
aid in source detection and size estimate, and subsequent con- 
strained multi-Gaussian fitting. This approach greatly increases 
the dynamical range between compact sources and diffuse emis- 
sion, irrespective of the local absolute value of the emission. The 
method is fully described elsewhere (Molinari et al., 2010a) and 
has been applied for this first attempt to generate source catalogs. 
As the thresholding for source detection is done on the curvature 
image (Molinari et al. , 2010a), the S/N of the detected sources is 
determined a posteriori measuring the ratio of the source peak 
flux over the rms of the residuals after the Gaussian fit. Source 
catalogs were generated for the two fields and for the 5 bands 
and are made available in tabular form in the online version of 
the paper. Catalogs completeness was estimated with artificial 
source injection experiments, and the peak flux levels for 80% 
completeness for the 70, 160, 250, 350 and 500/mi photometry 
are [0.5, 4.1, 4.1, 3.2, 2.5] Jy/beam for the /=30° field, and [0.06, 
0.9, 0.7, 0.7, 0.8] Jy/beam for the Z=59° field. The difference is 
entirely compatible with the very different intensity regimes of 
the underlying diffuse emission in the two fields. 

Estimating the source's physical properties requires that de- 
tection in the various band catalogs are merged in coherent 
SEDs, a process that can only be done coarsely in this early 
stage, but which is nonetheless useful for isolating 528 sources in 
the Z= 30° field and 444 sources in the /=59° field (see |Elia et al] 
|2010) l. The two observed fields encompass emission from re- 
gions at very different distances. In a considerable effort, which 
involved a critical re-evaluation of available data and evidence, 
and the collection of additional data for hundreds of previously 
unknown objects, Russeil et al. ( |2010 i provide recommended 
distances for a fraction of the detected sources (312 out of 
528, and 91 out of 444 sources for the two fields, respectively) 
for which the derivation of masses and luminosities is possi- 
ble. Adopting standard prescriptions for Class classification 
(LA>350iim/Lboi > 0.005, Andre et al.|20 00) results in almost the 
totality of sources being Class (90 out of 91 sources in /=59° 
and 306 out of 312 in /=30°, see |Elia et alpOTO] ). 



3.3. Filamentary star formation 

The most extraordinary feature exhibited by the Herschel maps 
is the ubiquitous pattern of filaments in the ISM structure. This 
is more apparent when we enhance the contrast of the filaments 
using the same method ( see §3.2|i as used for the source detection 



( |Molinari et al.| |2010aj). Here we start from the d 2 derivatives 
carried out in four directions (x, y, and the two diagonals), as for 
the standard detection method, and then create another image F 
of the same size so that the maximum curvature is selected for 
each pixel: Fjj = max[d 2 Xij,d 2 yij,d 2 D\jj,d 2 D2ij]. In this way 
we are following the direction of maximum curvature pixel-by- 




Fig. 3. Multidirectional second-derivative image (see text) of the 
Z=59° field at 250jum. The blue circles represent the compact 
sources detected at 250/mi. The filamentary structure of the ISM 
appears at various levels of intensity, i.e. curvature, and it is strik- 
ing how the detected compact objects are for the most part dis- 
tributed only along the brightest filaments. 



pixel for all compact features in the image. We show in figplthe 
result of this processing on the 1=59° field at 250/zm. 

The maps clearly show an interconnected maze of filaments 
at different levels of brightness (e.g. different levels of emis- 
sion intensity and curvature), and the striking aspect is that the 
compact sources detected at 250/mi are distributed for the most 
part along the brightest filaments. Interestingly, a similar sce- 
nario was also reported for Taurus by Goldsmith et al. (2008) 
where the physical conditions and spatial scales involved are 
radically different. Since the source integrated fluxes are esti- 
mated by fitting Gaussians on top of planar plateaus, the values 
of the local background at every wavelength are a by-product 
of our source extraction and, after applying the absolute cor- 
rection factors as recommended by Bernard et al. ( 2010[l and 
subtracting the foreground contribution estimated using Bohlin 



( 1975), can be used to estimate the local beam-averaged column 
density in the hosting filaments. The relationship between the 
mass of the detected cores (when the SED is reliable and the 
distance is known, see previous paragraph) and the local beam- 
averaged H2 filament column density is reported in Fig. HI The 
points for the /=59° field mostly lie in the range of column den- 



sities (ltr'cm^ < N(H 2 ) < 10 cm ) that corresponds to 
lmag < A v < lOmag, values that are entirely reasonable for 
the transition regime between diffuse ISM and dense molecular 
clouds ( |Cambresy|19"9"9"l|Snow & McC all 2006). Higher values 
of N(H2) are found for the points for the /=30° field, most likely 
due to the larger relative distances of the sources in this latter 
field. The core masses are spread between 1 and 10 4 M Q , with 
no indication of a correlation between the two quantities. The 
strong impression, however, is that of a threshold at Ay ~ lmag 
for the Z=59° field above which dense cores are found, a thresh- 
old that is evidently exceeded only in bright filaments (Fig. 13). 
More in particular, the Ay ~ Imag threshold corresponds to ~ 
17M pc~ 2 in molecular hydrogen , which is suprisingly close to 



the 10M Q pc value that Krumholz et al. (2009) find critical for 



the dust content in HI clouds to efficiently shield the cloud inte- 
rior from external FUV field and allow effective H2 formation. 
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Fig. 4. Mass of detected cores as a function of the hosting fila- 
ment local column density, for the cores where the distance is 
known and the mass could be estimated. The red diamonds rep- 
resent the cores in the /=59° region, while the plus signs repre- 
sent the cores from the /=30° region. The top X-axis represents 
the equivalent Ay. 



It is tempting to relate the appearance of clumps to an extinc- 
tion regime where the H1/H2 boundary shields the cloud interiors 
from interstellar FUV field, causing the photoelectric heating ef- 
ficiency to drop considerably and causing in turn a drop in dust 
and gas temperature ( |Tielens & Hollenbach 1985 1. This thresh- 
old value seems to be of the order ofAy~5-^10 for the 1=30° 
field most likely due to the relatively larger distances of sources 
in this field ( JRusseil et all|20T0"l l. 

The ubiquitousness of dense filaments in the ISM, the high 
degree of association between bright filaments and cores, and 
the suggestion of a column density threshold for the appearance 
of dense cores, all appear to coherently support a formation sce- 
nario that starts with the condensation of diffuse clouds into long 
filaments. As the column density increases, a threshold is ex- 
ceeded and denser star-forming (or potentially star forming) con- 
densations start to appear. A preliminary association with Spitzer 



24/im counterparts (Elia et al. 2010 1 suggests that our detected 



sources may be a mixture of protostellar and pre-stellar objects, 
although more work will be needed to ascertain the composition 
of this mixture. 

Since our source extraction also yields a measure of the core 
sizes, we are in the position to estimate their surface density 2. 
These values appear to be on average a factor 3-5 higher than 
the column densities of the underlying filaments as reported in 
Fig. m spanning a range of 0.03 g cirT 2 < £ < 0.5g cm" 2 , after 
changing units, with a mean value of 0.1 g cm" 2 . It is puzzling 
that the number of cores exceeding L/M ratios of a few (the ex- 
act number depending on the core masses), corresponding to the 
critical surface density threshold of 1 g cm -2 for the formation of 
massive stars ( |Krumholz & M cKee 2008), is not consistent with 



very few of the cores actually exceeding that critical threshold 
(see also Elia et al. 2010). The difference of a factor ~2 between 
the dust opacities that we used compared to Krumholz & McKee 
(2008JI is not sufficient to reconcile this apparent discrepancy. 
This result deserves more attention and needs to be confirmed in 
the future with more detailed and accurate analysis. 

Testing of large-scale "dynamical" star formation scenarios 
(e.g., Hartman n et al.|2001[ ), where filaments are formed in the 
post-shock regions of large Hi converging flows, is one of the 
original science goals of Hi-GAL. It is remarkable how the 
predictions from recent MHD numerical simulations ( jBanerjee] 
et al. 2009 \ of formation and subsequent fragmentation of fil- 
aments, agree with our results. Besides the morphological re- 
semblance of these simulations with the structures we see in our 
Herschel maps (Fig. IS), there is striking agreement of their pre- 
dictions with the N(H) regime we measured for our core-hosting 
filaments, as well as with the mass regime of the cores being 
formed. 

Instability and fragmentation of dense filaments has also 
been investigated in the context of helical magnetic fields en- 



closing the filaments by Fiege & Pudritz (2000); interestingly, 
the models predict the formation of regularly spaced condensa- 
tions at length scales that depend on the properties of the mag- 
netic field, the velocity dispersion, and density of the filament. 
The predicted length scale for filament velocity dispersion of 0.5 
km s" 1 and density of 10 4 cm" 3 is 2.8 pc, and curiously this is 
not at all far from the median distance of each source to its near- 
est neighbor: ~1.8pc for the 250//m sources in the /=59° field 
for the sources' average distance (Russeil et al. 2010). It is also 
interesting that the typical fragmentation length scale decreases 
with increasing filament density, in broad qualitative agreement 
with a higher spatial density of sources in the brightest filaments. 
This clearly deserves further investigation to be confirmed as a 
viable hypothesis. 



4. Conclusions 

The first science highlights presented in this paper, as well as 
in the accompanying papers in this volume and elsewhere, show 
that owing to its optimal use of unique Herschel characteristics 
of wavelength coverage, spatial resolution and mapping speed, 
the Hi-GAL survey has the potential to lead to a quantum leap 
in our understanding of large-scale Galactic star formation from 
cloud to cluster-forming clump formation and of the evolution 
of protoclusters and massive protostars. 

The outstanding feature emerging from these first images is 
the impressive and ubiquitous ISM filamentary nature. Dense 
cores seem to appear when a certain beam-averaged column 
density threshold is exceeded in close spatial association with 
these filaments. 
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